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First-Order  Perturbation  Corrections  to  the  Hartree-Fock  Approximation  for  Helium* 
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Variational  perturbation  theory  has  been  used  to  compute  the  first-order  wave  function  for  the  ground 
state  of  the  helium  atom,  in  a  scheme  where  the  Hartree-Fock  function  is  treated  as  the  xeroth-order  approxi¬ 
mation.  The  first-order  function  was  approximated  by  a  Hylleraas-ltke  expansion,  explicitly  containing  ru. 

We  obtain  a  total  energy  through  third  order  of  —2.90262  a.u.,  as  compared  with  the  nonrelativistic  eigen¬ 
value  of  —2.90372.  The  first-order  corrections  to  expectation  values  of  a  selection  of  one-  and  two-electron 
operators  have  been  computed.  The  calculations  have  also  been  done  for  the  iaoelectronic  ions,  H“,  Li+, 
and  Be**. 
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INTRODUCTION 

IT  is  well  known  that  the  Hartree-Fock  approximation 
neglects  the  details  of  the  interelectronic  repulsions, 
treating  them  in  only  an  average  fashion.  The  resulting 
error,  which  has  come  to  be  known  as  the  “correlation 
error,”  can  for  many  purposes  be  significant.  A  consider¬ 
able  body  of  literature  has  come  into  existence  dealing 
with  just  this  problem  in  a  variety  of  ways.1  Almost 
invariably,  the  approach  has  been  to  approximate  the 
many-electron  eigenfunction  by  a  variational  wave 
function,  which  explicitly  accounts,  to  some  extent  at 
least,  for  the  electron  correlation.  The  most  successful 
of  these  approaches  has  given  the  ground-state  energy 
of  the  helium  atom  to  well  within  the  experimental 
uncertainties.1 


02  In  this  paper,  we  follow  a  somewhat  different  ap- 
xdC  proach,  which  has  not  yet  received  much  computational 
attention,  although  the  formalism  has  been  discussed 
by  several  authors.*-*  The  total  Hamiltonian  is  sepa¬ 
rated  into  two  parts :  one  part  consisting  of  the  sum  of 
the  Fock  operators,  for  which  the  Hartree-Fock  orbitals 
are  eigenfunctions,  and  a  second  part  containing  every¬ 
thing  that  is  left  over,  which  is  treated  as  a  perturbation. 
The  iV-electron  Hartree-Fock  function  then  becomes 
the  first  term  (zeroth  order)  in  a  perturbation  expansion 
of  the  exact  wave  function,  all  the  correlation  effects 
being  thrown  into  the  higher  order  corrections. 

We  report  here  variational  calculations  of  the  first- 
order  function  for  the  ground  state  of  the  helium  atom 
(and  some  isoelectronic  ions),  based  on  the  variational 


principle  for  the  second-order  energy.7  In  particular, 
to  obtain  definitive  numerical  results,  the  first-order 
function  has  been  approximated  by  a  Hylleraas-like 
expansion  explicitly  involving  nonnegative  powers  of 
fit.  The  primary  purpose  of  this  investigation  is  to 
isolate,  for  a  simple  system,  how  much  of  the  correlation 
effects  are  contained  in  the  first-order  correction,  with 
the  intention  of  assessing  the  possibilities  of  the  scheme 
for  more  complex  systems. 

There  is  another  more  conventional  way  of  applying 
the  variation-perturbation  procedure  to  atomic  systems, 
which  deserves  mention  at  this  point.*  This  scheme, 
which  we  will  refer  to  as  the  "hydrogenic”  perturbation 
approach,  takes  the  sum  of  the  bare-nucleus  (hydrogen¬ 
like)  Hamiltonians  as  the  zeroth-order  Hamiltonian 
and  treats  the  entire  interelectronic  repulsion  as  the 
perturbation.  The  perturbation  expansion  of  the  energy 
and  all  expectation  values  becomes  a  descending  power 
series  in  the  nuclear  charge,  Z.  There  has  recently  been 
a  renewed  interest  in  this  approach,*-1*  and  we  will  be 
making  comparisons  between  it  and  the  present  “Har¬ 
tree-Fock  perturbation"  scheme. 

GENERAL  THEORY 

The  nonrelativistic  Hamiltonian  for  the  helium  atom 
is  given  by 

3C=A(l)-M(2)+l/r„,  (1) 

where  h(i)  is  the  “bare  nucleus”  Hamiltonian,1* 

h(i) - *A  i-Z/n.  (2) 


*  This  research  is  a  part  of  project  DEFENDER,  sponsored  by 
the  Advanced  Research  Projects  Agency,  Department  of  Defense, 
through  the  Office  of  Naval  Research. 

1  Rather  than  attempt  to  cite  all  the  relevant  literature,  the 
reader  is  referred  to  several  reviews  on  the  subject.  See,  e.g., 
P.  0.  Lowdin,  Advances  »'»  Chemical  Physics  (Interscience  Pub¬ 
lishers,  Inc.,  New  York,  1959),  Vol.  2;  J.  C.  Slater,  Quantum 
Theory  oj  Atomic  Structure  (McGraw-Hill  Book  Company,  Inc., 
New  York,  1960),  Vol.  2. 

*C.  L.  Pekeris,  Phys.  Rev.  112,  1649  (1958):  115, 1216  (1959). 

•C.  M(jller  and  M.  S.  Plesset,  Phys.  Rev.  46,  618  (1934). 

4  E.  M.  Corson,  Perturbation  Methods  in  the  Quantum  Mechanics 
oj  n-Electron  Systems  (Hafner  Publishing  Company,  Inc.,  New 
York,  1950),  pp.  150-156. 

•0.  Sinanoglu,  Proc.  Roy.  Soc.  (London)  A260, 379  (1961). 

4  M.  Cohen  and  A.  Dalgamo,  Pra^Phys.  Soc.  (London)  77, 
748  (1961) 

N0  OTi? 


For  the  ground  state,  the  Hartree-Fock  wave  function 
is  a  simple  orbital  product  function, 

_  *,=  *(1M2)  (3) 


1 H.  A.  Bethe  and  E.  E.  Salpeter,  Quantum  Mechanics  oj  One- 
and  Two-Electron  Atoms  (Academic  Press  Inc.,  New  York,  1957), 
pp.  122-123. 

'  E.  A.  Hylleraas,  Z.  Physik  65,  209  (1930). 

*  J.  Linderberg  and  H.  Shull,  J.  Mol.  Spectroscopy  5,  1  (1960). 

»  C.  Schwarts,  Phys.  Rev.  126,  1015  (1962). 

11 W.  B.  Somerville  and  A.  L.  Stewart,  Proc.  Phys.  Soc.  (London) 
80, 97  (1962). 

“  R.  E.  Knight  and  C.  W.  Scherr,  Phys.  Rev.  128,  2675  (1962). 

11  Atomic  units  are  used  throughout  this  paper.  Distance  is  in 
units  of  the  Bohr  radius,  and  energy  in  units  of  2Ru,  where  Ru 
is  the  reduced  mass  Rydberg. 
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where  we  have  suppressed  the  spin  part  of  the  function, 
and  the  orbitals  are  eigenfunctions  of  the  Fock  operator ; 


F(t)  v(t)=w(t), 

(4a) 

F(i)-h(i)+j,(i), 

(4b) 

fdV,(  Mj)|‘/r«). 

(4c) 

J,  is  the  usual  Coulomb  operator  of  self-consistent  field 
theory  and  i)  the  orbital  energy. 

The  Hamiltonian  (1)  can  be  rewritten  in  the  following 
fashion: 

3C-3C«+K  (S) 

with 

3e,-F(l)+F(2)}  (6a) 


F«l/r„-/,(l)-/,(2).  (6b) 


the  following  integral  an  extremum: 

/-<*>|(3Cs— «.)*,— 2(.,—  n*.>,  (13) 

with  respect  to  an  arbitrary  variation  of  leads  to 
Eq.  (9a).  For  the  exact  ¥i,  (13)  is  just  the  second-order 
energy,  while  for  an  approximate  function,  it  is  an  upper 
bound  to  «*  If,  now,  one  approximates  ¥i  by  a  linear 
combination  of  basis  functions, 

(14) 

and  applies  this  variational  principle,  the  matrix  ana¬ 
log  of  (9a)  results, 

(H|—  <«S)c»  «i«— v,  (IS) 

where  the  matrix  and  vector  elements  are  given  by 
(H»)*i**(^j>|UC#|$t),  (S)„-(*,|*t), 
(•),-<¥,|*>,  (▼),-(¥.!  F|*,>. 


Now,  the  Hartree-Fock  function  is  an  eigenfunction  of  For  the  present  calculations,  ¥,  was  approximated  by 
3C0|  which  we  treat  as  a  zeroth-order  Hamiltonian,  a  Hylleraas-like  expansion. 


3C^,-2#,,  (7) 

and  V,  the  “fluctuation  potential”,  becomes  a  perturba¬ 
tion.  The  exact  wave  function  and  energy  are  given 


by  the  perturbation  series, 

(8a) 

£=ft+«i+ft+ft-l - , 

(8b) 

where  the  ^<*s  are  solutions  of  the  equations 

(3Co—  ft)^i=  tiV'o—  Vit , 

(9a) 

(3Co—  tts)rh=  «*M-  V'h , 

(9b) 

and  the  energies  are  given  by  the  standard  expressions, 

«o=  2w , 

(10a) 

«i=<*.m*o> — (<p\j.\<p), 

(10b) 

ft=<*.m*i>, 

(10c) 

«»=(iM  v—  «i|fi). 

(10d) 

Also,  the  expectation  of  any  operator  /,  through  first 
order,  is  given  by 


¥1=  (8r)  "(piW+PiWW'. 

</» 

P-ff,  (17) 

where  f  is  an  additional  variational  parameter  available 
for  minimizing  h.  For  purposes  of  computing  the  matrix 
elements,  analytic  approximations  for  the  Hartree-Fock 
orbitals  were  used,  expanding  them  in  the  standard 
set  of  Slater-type  orbitals  (STO’s) ; 

«>-£  o«(2«)-+»(4w(2»)I)-‘'V-‘r-'.  (18) 

»,« 

All  the  matrix  elements  (16)  can  be  reduced  to  the  usual 
two-electron  integrals  of  the  conventional  variational 
calculations.  The  reduction  of  matrix  elements  involving 
J ,  to  these  integrals  is  given  in  an  Appendix. 

The  more  conventional  hydrogenic  perturbation 
scheme  is  obtained  by  taking  the  first  two  terms  of  (1), 
the  “bare  nucleus”  Hamiltonians,  as  the  zeroth-order 
Hamiltonian  and  the  entire  electronic  repulsion  as  the 
perturbation.  The  zeroth-order  solution  is  an  orbital 
product  of  hydrogenic  functions,  and  the  wave  function 
and  energy  expansions  become  descending  power  series 
in  Z,  the  energy  series  beginning  with  Z1. 


<*l/l*)i=<*>|/|*>>+2<*,|/|*,>.  (11) 

The  Hartree-Fock  total  energy,  in  terms  of  the  pertur¬ 
bation  energies,  is 

£hf=  ft+«i=2r(— (<e\Jr\ <p).  (12) 

It  should  be  noted,  in  Eqs.  (10)  and  (11),  that  we  have 
taken  ¥i  orthogonal  to  which  can  be  done  without 
loss  of  generality  and  insures  the  normalization  of  the 
total  wave  function  through  first  order. 

As  is  well  known,7  these  perturbation  equations  can 
be  obtained  by  a  variational  principle;  namely,  making 


RESULTS  AND  DISCUSSION 

The  energies  computed  for  helium  with  expansion 
lengths  as  high  as  50  terms  are  shown  in  Table  I.  The 
notation  Et  refers  to  the  total  energy  computed  through 
third  order,  i.e., 

£  «(• 

iH> 

The  second-  and  third-order  energies  appear  to  be  rea¬ 
sonably  well  converged  and  probably  represent  the  true 
values  to  at  least  four  decimal  places.  It  is,  however, 

■  •  ( 
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Table  I.  Hartree-Fock  perturbation  energiee  (In  a.u.) 
for  the  ground  state  of  helium. 


Maximum  No.  of 
of  i+j+u  terms 
InEq.(lf)  in*. 

H 

«« 

E, 

2 

7 

-0.04763 

0.00670 

-2.90262 

3 

13 

-0.04817 

0.00722 

-2.90263 

4 

22 

-0.04832 

0.00737 

-2.90263 

5 

34 

-0.04834 

0.00740 

-2.90262 

6 

50 

-0.04835 

0.00741 

-2.90262 

En-r- 

-2.86168 

An*- 

-2.90372* 

•See  Ret.  2. 


somewhat  discouraging  that  the  total  energy  is  no  closer 
to  the  eigenvalue.  The  error  is  about  0.03  eV,  which  may 
well  portend  substantial  energy  errors  if  the  scheme  were 
applied  to  much  larger  systems,  say  6  to  10  electrons. 

Our  numerical  experience  also  indicates  that  no  more 
than  the  number  of  figures  given  are  really  significant ; 
there  may  even  be  an  uncertainty  of  1  or  2  in  the  last 
place.  The  reason  is  apparently  the  approximate  nature 
of  the  zeroth-order  (Hartree-Fock)  functions,  which,  of 
course,  are  not  known  exactly.  This  is  a  limitation  that 
will  probably  apply  to  any  scheme  which  explicitly  in¬ 
corporates  the  Hartree-Fock  solutions  as  a  starting 
point.  It  is,  perhaps,  curious  to  note  that  the  total 
energy  through  third  order  is  quite  insensitive  to  the 
accuracy  with  which  is  approximated,  even  though 
the  individual  «i  and  «t  are  not. 

The  expectation  values  of  a  selection  of  operators 
and  their  first-order  corrections  are  shown  in  Table  II. 
The  first-order  corrections  were  computed  with  the  50- 
term  function;  calculations  with  the  smaller  wave 
functions  indicated  good  convergence  for  all  the  values 
given.  The  zeroth-order  part,  of  course,  is  just  the 
Hartree-Fock  expectation  value.  It  has  been  shown  that 
the  first-order  correction  to  a  Hartree-Fock  expectation 
value  of  any  one-electron  operator  should  vanish  identic¬ 
ally,'  and  this  is  borne  out  by  our  direct  calculations. 
This  feature,  in  fact,  provided  a  nice  additional  check 
on  our  numerical  work.  For  one-electron  operators, 
therefore,  the  Hartree-Fock  results  are  unchanged,  and 


Table  II.  First-order  corrections  to  the  Hartree-Fock  expecta¬ 
tion  values  of  some  one-  and  two-electron  operator!. 


/ 

(01/10) 

(01/1 1) 

(/>. 

</W 

rrl 

1.6874 

-0.0001 

1.6872 

1.6883 

ft 

0.9273 

0.0000 

0.9272 

0.9295 

fl* 

1.1849 

0.0000 

1.1849 

1.1935 

fir1 

1.0258 

-0.0483 

0.9291 

0.9458 

fit 

1.3622 

0.0351 

1.4325 

1.4221 

nf 

23698 

0.0859 

2.5416 

2.5164 

V  V» 

0.0 

-0.1060 

-0.2121 

-0.1591 

Tl'ft 

0.0 

-0.0429 

-0.0859 

-0.0647 

l(r,) 

1.7981 

-0.0002 

1.7977 

1.8104 

»('>•) 

0.1906 

-0.0574 

0.0759 

0.1063 

indeed  they  are  already  quite  good.  All  the  two-electron 
operator  expectation  values  are  substantially  improved, 
although,  just  as  with  the  energy,  the  corrections  may 
leave  something  to  be  desired.  The  delta  function, 
4(ri»),  for  instance  is  a  rather  severe  measure  of  correla¬ 
tion  effects,14  and  it  is  still  off  by  about  30%;  the 
Hartree-Fock  error  is  about  80%. 

The  vanishing  of  the  first-order  corrections  for  one- 
electron  operators  follows  from  a  general  “strong 
orthogonality”  property  of  the  first-order  function.*'4 
For  this  particular  separation  of  the  Hamiltonian,  the 
structure  of  the  first-order  Eq.  (9a)  is  such  that  ¥i  is 
automatically  insured  to  be  strongly  orthogonal  to  the 
orbital  <p,  i.e., 

/(I)-  j  dVtfx(  1,2)*(2)=0.  (19) 

This  integral  was  actually  computed  and  tabulated  for 
n  ranging  from  0.0  to  2.0  and  found  to  be  10~MO~* 
of  the  value  of  the  Hartree-Fock  orbital. 


Table  HI.  Energies  for  the  helium  isoeiectronic  sequence  (in 
s.u.) ;  comparison  of  the  hydrogenic  and  Hartree-Fock  perturbation 
results. 


H~ 

He 

Li+ 

Be»+ 

Ei  -0.48793 
Hartree-Fock  ^ 

E,  -0.51942 

-2.86168 

-0.04835 

0.00741 

-2.90262 

IIP 

NOON 

1  i  1 

-13.61130 

-0.04725 

0.00322 

-13.65533 

TTu-Hrnarntrl  ^1  “0.375 
Hydrogenic*  -0.52397 

-2.75 

-2.90332 

1  I 

3* 

—8 

-13.50 

-13.65549 

E^  -0.52775 

-2.90372 

-7.27991 

-13.65557 

First-order  corrections  were  also  computed  for  the 
isoeiectronic  ions  through  Z«=  4.  The  energy  results  are 
given  in  Table  III,  which  also  includes  the  analogous 
hydrogenic  perturbation  results.  It  should  be  noted  that 
the  50-term  first-order  function  for  H~  has  still  not  con¬ 
verged.  We  estimate  that  there  may  still  be  an  error  in 
the  energy  of  0.00015  a.u.,  which  however,  does  not  sub¬ 
stantially  alter  the  conclusion  that  the  first-order  func¬ 
tion,  in  an  isoeiectronic  sequence,  contains  an  increas¬ 
ingly  larger  portion  of  the  correlation  correction  as  Z 
increases. 

The  comparison  with  the  hydrogenic  perturbation 
results  are  somewhat  interesting.  Although  our  energy 
through  first  order  (Hartree-Fock  total  energy)  is 
better  than  the  corresponding  hydrogenic  energy,  the 
total  energy  through  third  order  is  worse.  While  the 
Hartree-Fock  approximation  provides  a  better  starting 
point,  the  hydrogenic  first-order  function  appears  to 
more  than  make  up  for  the  initial  deficiencies  of  the 
scheme,  energetically  at  least.  With  regard  to  expecta- 
tion  values  of  other  operators,  the  situation  generally 


•  See  Ret.  J. 


“  J.  Cooper  and  J.  B.  Martin,  Phyi.  Rev.  131,  1183  (1963). 
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appears  to  be  the  other  way  around.  Table  IV  lists  some 
expectation  values  computed  to  first  order  in  both 
schemes,  duplicating  some  of  the  numbers  of  Table  II. 
With  the  exception  of  rr1,  ru'1  (which  emphasize  those 
regions  of  space  important  for  the  energy),  and  t^u, 
the  present  results  are  closer  to  the  exact  values.  The 
results  for  one-electron  operators,  of  course,  merely 
reflect  the  superiority  of  the  Hartree-Fock  function. 

Some  calculations  were  also  done  on  helium,  using 
only  even  powers  of  ru  in  the  approximate  first-order 
function  (17).  This  is  equivalent  to  an  expansion  in 
Legendre  polynomials  of  the  angle  between  rt  and  u, 
i.e.,  it  is  stimulating  a  configuration  interaction  type 
approximation  to  the  first-order  function.  These  results 
are  given  in  Table  V,  where  the  5  and  S+P  notation 
refers  to  functions  with  no  rn  terms  and  functions  in¬ 
cluding  only  zero  and  second  power  terms,  respectively. 
While  we  were  not  able  to  use  sufficiently  large  expan¬ 
sions  to  definitely  determine  the  angular  limits,  the 
results  should  be  qualitatively  correct.  It  is  clear  that 


Table  IV.  Some  helium  expectation  values  through  first  order; 
comparison  of  the  hydrogenic  and  Hartree-Fock  perturbation 
calculations. 


Operator 

Hydrogenic* 

Hartree-Fock 

Exactb 

rr1 

1.687S 

1.6872 

1.6883 

u 

0.8203 

0.9272 

0.9295 

tf 

1.0488 

1.1848 

1.1935 

nr1 

0.9347 

0.9291 

0.9458 

fw 

1.3512 

1.4325 

1.4221 

f!^ 

2.0030 

2.5416 

2.5164 

Vi-V. 

-0.2666 

-0.2121 

-0.1591 

rrir 

-0.0473 

-0.0859 

-0.0647 

*(r.) 

1.6962 

1.7977 

1.8104 

0.0090 

0.0759 

0.1063 

most  of  the  first-order  correction  is  due  to  angular 
effects. 

In  all  these  calculations,  a  five-term  representation 
of  the  Hartree-Fock  orbital  was  used.  In  particular,  it 
was  expanded  in  five  Slater-type  lr  orbitals  (fully 
optimized).  To  determine  the  effect  of  using  a  cruder 
approximation  to  the  Hartree-Fock,  the  calculations 
were  repeated  with  a  three-term,  “single  exponential” 
expansion  of  the  H-F  orbital.  The  results  for  <i  and  <s 
differed  by  4  and  5  units,  respectively,  in  the  fifth  place 
from  our  best  values  in  Table  I. 

CONCLUSIONS 

The  present  pilot  calculations  on  two-electron  atoms 
suggest  that  a  perturbation  scheme  based  explicitly  on 
the  Hartree-Fock  model  will,  in  first  order,  account  for 
a  major  share  of  the  correlation  effects.  The  energy  and 
expectation  values  (with  the  exception  of  one-electron 
properties)  are  substantially  improved  by  the  first-order 


Table  V.  Angular  limit  calculations  for  helium. 


Angular 

type 

No.  of 
terms 
in*i 

•t 

«s 

E, 

S 

16 

-0.01794 

0.00083 

-2.87879 

S+P 

28 

-0.04253 

0.00413 

-2.90007 

All  types 

50 

-0.04835 

0.00741 

-2.90262 

correction.  Hence,  for  atomic  and  molecular  properties, 
this  procedure  may  provide  a  reasonably  good  “second 
stage  of  approximation”  after  the  Hartree-Fock. 

However,  if  one  is  interested  in  a  high  degree  of 
absolute  precision,  some  caution  may  be  in  order.  This 
procedure  obtains  97.5%  of  the  correlation  energy  for 
helium,  corresponding  to  an  absolute  error,  due  to 
higher  order  effects,  of  0.03  eV.  Since  this  error  will 
almost  certainly  increase  with  the  number  of  electrons, 
a  first-order  treatment  may  well  become  inadequate  for 
much  larger  systems  than  helium.  Our  E~  results  also 
indicate  that  the  higher  order  effects  of  correlation  are 
most  important  for  negative  ions,  and  they  would  have 
to  be  explicitly  taken  into  account  to  get  reliable  results 
for  such  properties  as  electron  affinities. 

In  order  to  obtain  definitive  quantitative  results,  we 
have  explicitly  introduced  the  interelectronic  coordinate, 
rn,  into  the  first-order  trial  functions.  In  attempting  to 
extend  this  technique  to  larger  systems,  however,  one 
gets  the  same  three-electron  integral  problems  as  in 
the  usual  variational  approaches;  they  are  brought  in 
by  the  exchange  effects  in  the  Fock  operator.  To  avoid 
these  integral  difficulties,  one  would  presumably  have  to 
make  use  of  some  kind  of  configuration  interaction  type 
approximation  for  *i,  and  it  may  be  of  some  interest 
to  examine  the  convergence  properties  of  this  method 
when  applied  to  this  perturbation  scheme.  Its  generally 
slow  convergence,  however,  seems  to  warrant  a  rather 
pessimistic  outlook.u>u 
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APPENDIX 

The  matrix  elements  of  Eq.  (16)  are  readily  reduced 
to  sums  over  the  conventional  two-electron  integrals 
that  always  arise  in  any  variational  calculation  on 
helium,  which  uses  rn  explicitly.  The  only  operator 
that  warrants  any  discussion  is  the  Coulomb  operator 
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